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Abstract. We present algorithms used in the computational part of the article "Special homogeneous 
linear systems on Hirzebruch surfaces" . 



The main aim of this paper is to give the detailed description of algorithms used in proving the main 
Theorem 6 in |Dum 09a| . We will not repeat the definitions presented in |Dum 09a] . 

All algorithms are presented in self-explaining pseudo-code. We will use indentation to make our 
algorithms easier to read. The < — means an assignment, i.e. in the line 

A< — B; 

we force A to be equal to B. We use only integers, so m > 2 means that m G {2, 3, 4, . . . }. The command 
return finishes our algorithm immediately, so further commands (if any) won't be executed. The control 
structures (if, for, repeat and so on) will be used in two versions. The first one, with only one command 
executed: 

if . . . then command; 

and the second one, with possibility of more that one command to execute: 
if . . . then 

command A; 

command B; 

. . . , 
end if 

1. Basic algorithms 

The first algorithm is used to m-reduce a given diagram, according to [Dum 09al Definition 20] . 
Algorithm reduce 

Input: TO > 2, a diagram D = c)ia0(ai, . . . , a^). 
Output: c)iag(6i, . . . ,bk) — Tedm{D) or not reducible if D is not m-reducible. 

if fc < TO then return not reducible; 
for j = 1, . . . , k do bj < — aj ; 

U < — 0; 

for j — k, . . . ,k — m + 1 do 

if Gj < TO then r < — aj else r < — max({l, . . . , to} \ U); 



aj - r- 



U < UU{r}; 

end for 

if J7 = {1, . . . , to} then return 5iag(6i, . . . ,bk) else return NOT REDUCIBLE; 



Example 1. Let us compute REDUCe(3, c)ia0(5, 5, 4, 2)). First we put (&i, 62, ^3, bi) = (5, 5, 4, 2) and take 
U = 0. Next, for j = 4, 3, 2 we will proceed as in the "for" loop. For j = 4 we have 04 = 2 < 3 = to, 
so we take r = 2, 64 = 2 — 2 = 0, and we add {2} to the set U. The second step, for j — 3, gives 
r = max({l, 2, 3} \ {2}) = 3, so r = 3 and b^ — 1. Before passing to the third step, we put U = {2, 3}. 
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In the last step we will have r ~ max({l, 2, 3} \ {2, 3}) = 1, so 62 = 4. At the end we have U = {1, 2, 3}, 

so reduce(3, J)iafl(5, 5, 4, 2)) = Oiag(5, 4, 1). 

We will apply a sequence of reductions to a diagram D, so we define an auxiliary algorithm SEQUENCE- 
REDUCE. By red^^ we denote 

red^^ = redmo-- -ored^. 

k 

Algorithm SEQUENCE-REDUCE 

Input: m > 2, a number A; > 1, a diagram D. 
Output: redl^\D) or NOT REDUCIBLE if reduction fails at some step. 

repeat k times 

D i — REDUCE(m, D); 

if D = NOT REDUCIBLE then return NOT REDUCIBLE; 
end repeat 
return D; 

Example 2. Let us compute sequence-reduce(4, 3,i)ia£|([4]^^, [5]^^)). The reductions goes as follows: 
i)iafl([4] 5, 5, 5, 5, 5) i)iafl([4] 5, 4, 3, 2, 1) OtO0(4, 4, 4, 5) i)iafl(3, 2,1,1), 

so 

SEQUENCE-REDUCE(4,3,3iafl([4]^^, [5]^^)) = 0iafl(3,2, 1, 1), 

while 

SEQUENCE-REDUCE(4,4,0tOg([4]^^, [5]^^)) = NOT REDUCIBLE. 

The following two algorithms will be used in the algorithm CH. The first one simply reduces all 

diagrams from a given set (ommiting not reducible ones) , the second finds all diagrams from a given set 
V which cannot be reduced to some diagram from the second given set Q. 

Algorithm RED 

Input: m > 2, a number fc > 1, a set I? of diagrams. 
Output: the set Q = {veS^^ [D) : D £T>,D is m-reducible k times}. 

0^0; 

for each D € V do 

D < SEQUENCE- REDUCE(rn, k,D); 

if D ^ NOT REDUCIBLE then G < — G U {D}; 
end for each 
return G', 

Algorithm REDOUT 

Input: m > 2, a number fc > 1, a set P of diagrams, a set G of diagrams. 
Output: the set {D e D : red^'(D) ^ 5}. 

for each D do 

G < SEQUENCE-REDUCE(m, k, D); 

itG eG then V< — D \ {D}; 
end for each 
return V; 

We will often reduce a diagram D as many times as possible, so we define an auxiliary algorithm 

TOP-REDUCE. 

Algorithm TOP-REDUCE 

Input: TO > 2, a diagram D ~ diaQ{ai, . . . , at). 
Output: G = Diafl(6i, . . . ,bk) = redm(redm(- • • (-D) . . . )) such that G is 
not TO-reducible. 
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repeat 

G < — REDUCE(m, £)); 
if G = NOT REDUCIBLE then return D; 
D< — G; 
end repeat 



Observe that if D is not m-reducible then top-REDUCe(£') = D. 

Example 3. Let us compute top-REDUCe(3, J)ia0(5, 5, 4, 2)). We have 

reduce(3, diag{5, 5, 4, 2)) = c)ia0(5, 4, 1), 

reduce(3, c)ia0(5, 4, 1)) = 5ia0(3, 1) 

reduce(3, I)iag(3, 1)) = not reducible. 

Hence top-reduce(3, Oia0(5, 5, 4, 2)) = t)ia0(3, 1). 

Now we present an algorithm to find all admissible /i-c)iag(6i, . . . , 6m_i)-tails for multiplicity m, see 
jPum 09al Definition 35]. By the length of a diagram D — 5iag(ai, . . . , Ofe) we denote the number of its 
non-zero layers, 

length(i?) = #{j : a, > 0}. 

Algorithm h-tails 

Input: m > 2, h > m, D = (3\ag{bi, . . . , &m-i). 
Output: the set V of all admissible h-D-ta.iis for multiplicity m, 
or ERROR if some reduction stops too early. 

V < — 0; 
repeat 

V< — VU{D}; 

D < — diagih) + D; 

D < TOP-REDUCE(m, D); 

if length(D) > m then return error; 

\i D eT) then return V; 
end repeat 



Example 4. We will find all admissible 4-c)iag(0, 0)-tails for multiplicity 3. The set V is empty at the 
beginning, and we put D — i)iag(0, 0) ~ into V. Now we take new D — c)ia0(4), reduce it as many 
times as possible, but, in fact, top-REDUCe(3, c)iag(4)) = c)iag(4). So we go at the beginning of the 
"repeat" loop and add c)ia0(4) to the set V. Now we take new D = l3\aQ{4) + t)ia0(4) = c)iog(4, 4), still 
it cannot be 3-reduced. Thus we add it into V, which now is equal to {0, 5iag(4), Oia0(4, 4)}. Taking 
D = c)io0(4) + c)ia0(4, 4) = 5iag(4,4, 4) we obtain top-REDUCe(3, t)iag(4, 4, 4)) = 0. Since the last 
diagram already belongs to 2?, the algorithm terminates. 

Observe that the size of V depends also on h and t)iag(5i, . . . , 5m_i). For example, 

h-tails(3, 5, c)ia0(0, 0)) = {0, c)ia0(5), c)ia0(5, 5), c)ia0(3), I)ia0(5, 3), 

0100(4, 3), c)ia0(4, 2), t)ia0(4, 1), Oia0(3, 1)}, 

while 

#H-TAILS(3, 5, ()iO0(l, 0)) = #H-TAILS(3, 5, 0100(0, 0)) + 4. 

Example 5. We will compute h-tails(5, 9, J)ia0(O, 0, 0, 0)). The computations can be written in the 
following short form: 

— > Oia0(9) — > t)ia0(9, 9) — > c)ia0(9, 9, 9) — > c)ia0(9, 9, 9, 9) — > 

T)iag{9, 9, 9, 9, 9) — > l>iag{8, 7, 6, 5, 4) — > diag{7, 5, 3) — > c)ia0(9, 7, 5, 3) — > 

J)ia0(9, 9, 7, 5, 3) — * diagiS, 7, 3) — > Uag{9, 8, 7, 3) — » diag{9, 9, 8, 7, 3) — > 

Mo0(8, 7, 4, 2) — > c)ia0(9, 8, 7, 4, 2) — > diag{8, 5, 2) — > diag{9, 8, 5, 2) — > 

Oio0(9, 9, 8, 5, 2) — > c)ia0(8, 6, 4) — > i)ia0(9, 8, 6, 4) — > c)ia0(9, 9, 8, 6, 4) — > 

5ia0(8, 7, 5, 1) — > c)ia0(9, 8, 7, 5, 1) — > diag{7, 5, 3) 



3 



and wc finish, since the last diagram has been found earlier. 



Since we are interested in collecting all admissible h-D-tai\s for all D G V, we present an auxiliary 
algorithm, called ltails. 

Algorithm ltails 

Input: TO > 2, /i > m, a set P of diagrams. 
Output: the set Q of all admissible ft,-Z?-tails for multiplicity to and D G T>, 
or ERROR if some reduction stops too early. 

for each D £ T> do 

H < — H-TAlLS(m, h, D); 

iiTi. — ERROR then return error; 

end for each 
return Q; 

Sometimes we also want to find all (top)reductions of diagrams of the form 

diasiih]''") + D : D e V 
for some fixed n and a set T) of diagrams. 
Algorithm atails 

Input: m > 2, h > m, n > 0, a set V of diagrams. 
Output: the set G of all diagram from {Uag{[h]^") + D : D e V} 
reduced as many times as possible, 
or ERROR if some reduction stops too early. 

g^0; 

for each D E T> do 

G^Oiag([/i]X") + Z); 

G < TOP-REDUCE(m, G); 

if length(G) > to then return error; 

end for each 
return G; 



Example 6. Let us compute atails (3, 5, 2, 2?) for 

V ^ {0,5ia0(5),c)ia0(5,5),t)ia0(3)}. 
For each D gV we execute top-REDUCe(3, Oiag(5, 5) + £>). The result is 

{t)ia0(5, 5), t)ia0(3), 5ia0(5, 3), c)ia0(4, 3)}. 
Observe that ATAILS(3, 4, 2, 2?) — error, since 

top-reduce(3, Oia0(4, 4, 5)) = c)ia0(4, 2, 2), 

which is too long. 

Now our aim is to enumerate all admissible c)ia0(ai, . . . ,am)-tails for multiplicity to, see 'Dum 09a| 
Definition 39]. All admissible tails could be found by iterating symbolic reductions, see Dum 09al 
Definition 37] and the discussion after Example 40. Therefore we present an algorithm to produce all 
symbolic reductions of a given diagram of the form Oiag(ai, . . . ,am, [x]^*^). This amounts to substitute 
jj^jxfc g^Y[ reasonable integers and reduce obtained diagrams. Observe that for A; = 0, the symbolic 
reduction is equal to the TO-reduction. 

If D = t)ia0(ai, . . . , ak) then by cut(Z?, r) we denote the diagram given by 

A hia2{ai,...,ar) r<k, 
^'^'^''^'^ = [0 r>k. 



Algorithm SYMB-REDUCE 



Input: m > 2, i)ia£|(ai, . . . ,am, [x]^*^), < fc < m — 1. 
Output: the set V of all symbolic reductions of Ota0(ai, . . . , am, [x]^*^). 

Vi — 0; 

if A; = then 

D < — REDUCE(m,OtO0(ai, . . . ,0™)); 

if D ^ NOT REDUCIBLE then V < — {£>}; 

return V; 
end if 

for each (ci, . . . , Ck) satisfying min{m + 1, am} > ci > C2 > • • • > Cfe do 
D < — £iia0(ai, ... ,0^,01, . . . ,Cfc); 
D < — REDUCE(m, £)); 
if D 7^ NOT REDUCIBLE then 
G < — cut(£', m); 
e < — lengtli(£)) - length(G); 
G + c)io0([x]x^); 
— VLI{G}; 
end if 
end for each 
return V; 

Example 7. Let us compute symb-reduce(3, Diag(5, 5, 5, x, x)). In the "foreach" loop we must consider 

all pairs (01,02) satisfying 4 > ci > 02- For each such pair wc take Oiag(5, 5, 5, ci, C2), 3-rcduce it, and 
change the fourth and fifth number into symbols x, if necessary. We present the computations in the 
following table. 



(C1,C2) 


t)ia0(5,5,5,ci,C2) 


red3(OtO0(5,5,5,ci,C2)) 


the result 


(4,4) 


0100(5,5,5,4,4) 


Oia0(5,5,4.2, 1) 


Oio0(5, 5, 4, X, x) 


(4,3) 


Mas(5,5,5,4,3) 


Oia0(5,5,4,2) 


Oio0(5,5,4,x) 


(4,2) 


t)ia0(5,5,5,4,2) 


()ia£i(5,5,4,l) 


0x00(5, 5,4, x) 


(4,1) 


0100(5,5,5,4,1) 


Oio0(5,5,3, 1) 


Oio0(5, 5, 3, x) 


(4,0) 


£)ia0(5,5,5,4) 


Oia0(5,4,3,l) 


Oio0(5,4, 3,x) 


(3,3) 


t)ia0(5,5,5,3,3) 


Oia0(5,5,4,l) 


0x00(5, 5,4, x) 


(3,2) 


0100(5,5,5,3,2) 


Oia0(5,5,4) 


0x00(5,5,4) 


(3,1) 


t)ia0(5,5,5,3,l) 


Oio0(5,5, 3) 


Oio0(5,5, 3) 


(3,0) 


0iag(5, 5,5,3) 


0X00(5,4,3) 


0X00(5,4,3) 


(2,2) 


0100(5,5,5,2,2) 


NOT REDUCIBLE 




(2,1) 


Oio0(5,5,5,2,l) 


OiO0(5,5,2) 


Oio0(5,5,2) 


(2,0) 


Oio0(5,5,5,2) 


0X00(5,4,2) 


0X00(5,4,2) 


(1,1) 


Oio0(5, 5, 5, 1, 1) 


NOT REDUCIBLE 




(1,0) 


Oia0(5,5, 5,1) 


OiO0(5,3,2) 


Oio0(5,3,2) 


(0,0) 


Oia0(5,5,5) 


0X00(4, 3, 2) 


0X00(4, 3, 2) 



Now we can enumerate all admissible Oio0(ai, . . . ,am)-tails for multiplicity m. To do this, we will 
consider all symbolic reductions of symbolic reductions of . . . and so on. Each diagram obtained in this 
way, which is short enough (i.e. without x and with length at most m—1), satisfies the desired property. 

Algorithm TAILS 

Input: m > 2, Oio0(ai, . . . , ttm). 
Output: the set V of all admissible Oxo0(ai, . . . , am)-tails. 

Vi — 0; 

{Oxo0(ai,...,a™,[x]x(— 1))}; 

repeat 

choose W gW; 

W < — W\{W}; 

TZ < — symb-reduce(to, W); 
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W < — WUTl; 

for each D € TZ do 

if length(L») < m then V < — VU {D}; 
end for each 
until W = 



To show that the above algorithm terminates after a finite number of steps, observe that in each step, 
after choosing W eW and producing TZ — SYMB-REDUCE(m, W), we have 

#iy > max{#i? : Ren}, 

where 

#I)ia0(ai, . . . ,afc, [x]*^) = ai H h Ofe. 

Example 8. We will find tails(2, Oia0(5, 5)). The example for m — 3 would be a bit too long. The idea 
is to compute consecutive symbolic reductions of c)iag(5, 5, x). In the first step we proceed as in Example 

El 



D 


Cl 


Die) 


red2(i?(ci)) 


the result 


5ia0(5,5,x) 


3 


c)ia0(5, 5,3) 


c)iag(5,4,l) 


J)ia0(5, 4, x) 


Oia0(5,5,x) 


2 


c)ia0(5,5,2) 


fia0(5,4) 


fia0(5,4) 


()ta0(5,5,x) 


1 


c)ia0(5, 5, 1) 


^ia0(5,3) 


i)ia0(5,3) 


5ia0(5,5,x) 





diag{5, 5) 


fia0(4,3) 


fia0(4,3) 



In the next step we take all obtained diagrams and perform all possible symbolic reductions. 



Again, in the third 



D 


Cl 


Die) 


Yed2{D{ci)) 


the result 


i)ia0(5,4,x) 


3 


c)ia0(5,4,3) 


Olo0(5,3,l) 


Oio0(5,3,x) 


5io0(5,4,x) 


2 


c)ia0(5,4,2) 


0100(5,3) 


Oio0(5,3) 


5ia0(5,4,x) 


1 


c)ia0(5,4,l) 


Oio0(5,2) 


Olo0(5,2) 


t)ia0(5,4,x) 





5100(5,4) 


Oio0(4,2) 


Olo0(4,2) 


fia0(5,4) 




5100(5,4) 


Oio0(4,2) 


Olo0(4,2) 


5ia0(5, 3) 




c)la0(5,3) 


Oio0(4, 1) 


Olo0(4, 1) 


5io0(4, 3) 




0100(4,3) 


0100(3, 1) 


Oio0(3, 1) 


ep: 

D 


Cl 


Die) 


red2(i?(ci)) 


the result 


5ia0(5,3,x) 


3 


0100(5,3,3) 


Olo0(5,2,l) 


Oio0(5,2,x) 


5ta0(5,3,x) 


2 


0100(5,3,2) 


Olo0(5,2) 


Olo0(5,2) 


()ia0(5,3,x) 


1 


0100(5,3,1) 


Olo0(5, 1) 


Oio0(5, 1) 


£)ia0(5,3,x) 





Olo0(5,3) 


Olo0(4, 1) 


Olo0(4, 1) 


5io0(5,3) 




Oio0(5,3) 


0100(4, 1) 


Olo0(4, 1) 


fia0(5,2) 




Olo0(5, 2) 


Olo0(4) 


Olo0(4) 


5ia0(4,2) 




Olo0(4, 2) 


Olo0(3) 


Olo0(3) 


£)ia0(4, 1) 




Olo0(4, 1) 


Olo0(2) 


Olo0(2) 


5ia0(3, 1) 




Olo0(3, 1) 


Olo0(l) 


Olo0(l) 



The diagrams with length 1 are no more 2-reducible, so they won't produce any additional admissible 
tail. We present the fourth step: 



D 


Cl 


Die) 


red2(Z?(ci)) 


the result 


Olo0(5, 2,x) 


2 


Olo0(5,2,2) 


Olo0(5, 1) 


Oio0(5, 1) 


Olo0(5,2,x) 


1 


Olo0(5,2,l) 


Olo0(5) 


Olo0(5) 


Oio0(5, 2,x) 





Olo0(5,2) 


Olo0(4) 


Olo0(4) 


Olo0(5,2) 




Olo0(5,2) 


Olo0(4) 


Olo0(4) 


Olo0(5, 1) 




Olo0(5,l) 


Olo0(3) 


Olo0(3) 


Olo0(4, 1) 




Olo0(4, 1) 


Olo0(2) 


Olo0(2) 



In the final step we must reduce Olo0(5, 1) to obtain Olo0(3). We collect all obtained diagrams of length 
at most 1, thus 

TAiLS(2,Oio0(5,5)) = {Olo0(5),Olo0(4),Olo0(3),Olo0(2),Oio0(l)}. 
Observe that we can skip some of the above reducing. 
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2. Algorithms to compute sets V 



In [Dum 09al Section 7] we construct various sets of diagrams. Each set serves for showing that some 
given family of systems contains only non-special ones. To be more precise, we define a family S of 
systems together with a finite set V of diagrams such that if for each D G T) and r — [t^^ttJ the 

( 2 j 

systems C{D;m^'^) and £(_D; m^'-'"+^-') are non-special then S contains only non-special systems. 

The first algorithm computes the set V from [Dum 09al Proposition 42]. For D — J)iag(ai, . . . , ak) let 
rev(_D) = Oia0(afe, . . . , ai). Similarly, for a set V of diagrams, let 

rev(X») = {Tev{D) : D eV}. 

Algorithm setbign 

Input: m > A, N > m. 
Output: the set V from Proposition 42. 

V < 0; 

C < — H-TAiLS(m,TO+ l,Oia0([O]"=("'-i))); 
for j ~ m + 2, . . . , 2m — 3 do 

C < — ATAiLS(m, j, N, £); 

£ < — LTAiLS(m, 
end for 

C < — ltails(to, 2m - 2, C); 

n < — TAiLS(TO,c)ia0([2m - 1]''")); 

for each {L, R) e C x TZ do 

V< — VU {rev(L) + t)iafl([2TO - 2] + R}; 
end for each 
return V; 

Example 9. We will show the example for m = 5, = 11 (which is a part of our computation to prove 

Theorem 6 in |Dum 09a] V We will not present all the details, since the output would be too big. In our 

case we do the following: 

C < — H-TAILS(5, 6, 0100(0, 0, 0, 0)); 

C < — ATAILS(5, 7, 11, £); 

C < — LTAILS(5,7, £); 

C < — LTAILS(5, 8, £); 

n i — TAILS(5, Viag{9, 9, 9, 9, 9)); 

In the first step we obtain 

C = {0, c)iag(6), c)ia0(6, 6), c)iag(6, 6, 6), c)ia0(6, 6, 6, 6)}. 

In the second step, for each D e C, wc take top-REDUCe(5, t)iag([7] ^^^) + D). After reducing, we will 
have 

C = {c)ia0(7, 6, 4), c)ia0(7, 7, 6, 3), Oiafl(6, 4, 3, 1), Oia0(5), t)ia0(7, 4)}. 
In the third step we look for top-REDUCE of all diagrams of the form 

5ia0( [7]^ '=)-!-£', k>0,DEC. 

We will not enumerate all of them, since after this step, = 53. In the fourth step we look for 

J)iafl([8]'"') ^-i:*, k>0,DeC, 

and the resulting set contains 119 diagrams. Now we look for admissible tails. After computations, we 
obtain TZ of cardinality 147. 

Now we must "glue" diagrams from C and TZ to produce 119 • 147 = 17493 diagrams in V. One can 
check that, for example, 

Uag{8, 6, 3, 1) e C, <)iag{7, 6, 5, 4) e TZ, 

so we have 

0100(1,3,6,8, [8]''",7,6,5,4) e P. 
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Our next algorithm computes the set T) from Proposition 44 in [Dum 09a| . 
Algorithm setbign23 

Input: 2 < m < 3, iV > m. 
Output: the set V from Proposition 44. 

V < — 0; 

C < — H-TAiLS(m,m+ l,i)iafl([0]^("'-i))); 

C < — ltails(to,to + 2,C); 

TZ< — TAiLS(TO,c)ia0([m + 3]^™)); 

for each {L, R) £ C x TZ do 

V< — VU {rev(L) + diaQ{[m + 2]'"^) + R}; 
end for each 
return V; 

Example 10. The example for m — 3 would be very nice and illustrating, but also too long. So we will 

compute SETBIGn23(2, 2). We have three steps: 

C < — H-TAiLs(2,3,c)ia0(O)); 

C < — LTAILS(2,4,/:); 

n < — TAILS(2, f)ia0(5, 5)); 

In the first step we obtain 

£ = {0,c)iag(3)}. 

In the next step 

£ = {5100(4), ma0(3),Oia0(2),c)ia0(l)}. 
In Example [8] we have shown that 

TZ = {c)ia0(5), t)ia0(4), c)ia0(3), Oia0(2), c)ia0(l)}. 

So the final set is 

V = {Oia0(a, [4]^^6) : 1 < a < 4, 1 < 6 < 5}. 
Our next algorithm computes the set V from [Dum 09al Proposition 47]. 
Algorithm setbignb 

Input: m > 2, N > m, b > m + 2. 
Output: the set V from Proposition 47. 

V < — 0; 

C < — h-tails(to,to+ l,Oia0([O]^(™-i))); 
for j = TO + 2, . . . , 6 — 1 do 

C < — ATAILS(to, j, N, £); 

C < — ltails(to, j, £); 
end for 

C < — ltails(to, b, C); 

n< — h-tails(to,6+ l,ma0([O]^('"-i))); 

for each {L,R) e C x TZ do 

V< — VU {rev(L) + t)ia0([6] + i?}; 
end for each 
return V; 

Example 11. We will show the example for m = 5, N = 11, b ~ 8 (which is a part of our computation 

to prove Theorem 6 in [Dum 09aj ). In our case we do the following: 

C < — H-TAILS(5, 6, diagiO, 0, 0, 0)); 

C < — ATAILS(5, 7, 11, £); 

C < — LTAILS(5, 7,£); 

C < — LTAILS(5,8,£); 
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n i — H-TAILS(5, 9, 5ia0(O, 0, 0, 0)); 



This is very similar to what we did in Example[21 except for the last step. So we will have = 119. 
Now we must look for all possible top-reductions of t)iag([9] ^'^). By Example [5] we have ^TZ = 15, so 
#2? = 15-119 = 1785. 

The next algorithm computes the set T) from [Dum 09al Proposition 50]. For a diagram D = 
c)ia0(ai, ...,ak) let 

\diaQ{ak-e+i,ak-e+2, ■ ■ ■ ,ak) ^ < k, 
\d i>k. 

Algorithm SETNB 



cut 



Input: m > 2, n > 2, B > 2m - 1. 
Output: the set V from Proposition 50. 

V < — 0; 

G ^ mafl([m + 1] , [B]'''\B + 1); 

H < — cut(G', m); 

K < — cut(G, n{B ~m) -m + 1); 

TZ < — TAiLS(m, H); 

for each R (£ TZ do 

V < — VU{K + R}; 
end for each 
return V; 



Example 12. Let us compute SETNb(3, 2, 6). We wih have 

G = Oia0(4,4,5,5,6,6,7), 

so we take 

iJ = 5ia0(6,6, 7), = Dia0(4, 4, 5, 5). 

Now we must find TZ — TAiLS(m, H) by considering all symbolic reductions of t)iag(6, 6, 7, x, x), and take 

V = {Oia0(4, 4,5,5) + R: Re TZ}. 

The next algorithm computes the set V from [Dum 09al Proposition 52]. 
Algorithm setnba 

Input: m > 2, n > 2, b > m + 1, A > 0. 
Output: the set V from Proposition 52. 

V < — 0; 

TZ < — H-TAiLS(m,5+ l,Oiafl([0]"=("'-i))); 
for each R E TZ do 

V< — 2? U {c)iafl([TO + 1]>^", . . . , [&]><", [b + l]x('4+i)) + R}- 
end for each 
return T>; 



Example 13. For (m, n, b, A) — (3, 2, 5, 0) we will have 

TZ = H-TAILS(3, 6, 5iafl(0, 0)) 

= {0, t)io0(6), 5ia0(6, 6), c)ia0(4, 2), Oia0(5, 1)}, 

so 

D = {Oia0(4, 4, 5,5,6) + R: Re TZ}. 
The next algorithm computes the set T> from [Dum 09al Proposition 54]. 
Algorithm SETPB 
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Input: m>2,B> 3{m - 1). 
Output: the set V from Proposition 54. 



V < — 0; 

7^ < — TAiLS(r7i, DiagiB - m + 2, B - m + 3, . . . , B + 1)); 
for each i? G 7?. do 

V< — VU {t)iafl(l, 2, . . . , B - m + 1) + i?}; 
end for each 
return V; 



Example 14. We will compute SETPB(3, 9). Hence our computation starts with 

n = TAiLs(3,I)iag(8,9, 10)). 
We obtain ^TZ = 28 and take 28 diagrams of the form 

mafl(l,2, ...J)+R.: Ren. 
For example, we wih have c)ia0(l, 2, 3, 4, 5, 6, 7, 6, 5) G V. 

The next algorithm computes the set T) from I Dum 09al Proposition 56]. 
Algorithm setpba 

Input: m > 2, B > m, A > b. 
Output: the set V from Proposition 56. 

V < — 0; 

n < — H-TAiLS(m,6+ l,^iafl([0]^('"-i))); 
for each R € TZ do 

V^VU {dias{[b + l]x('4+i)) + R}- 
end for each 
return V; 



Example 15. It is easy to check that 

SETPBA(3,7,7) = {5iQ0([8]><«),c)ia0([8]><9)^j,i^g([8]xiO)^j,i^g([8]x8^5^1)^ 

c)ia0([8] 6, 2), c)ia0([8] 7, 3), 5ia0([8] 7, 5)}. 

3. Checking non-speciality 

We begin with an algorithm to decide whether a given system C{D;m^^) is special or not. The 
computations will be performed over ¥p and for some randomly chosen coordinates of points. Therefore, 
if our "specialized" system is non-special then obviously the general one is also non-special. In the 
opposite case we only know that our method does not work. We begin with preparing the matrix for our 
system (see |Dum 09c] ). 

Algorithm i-MATRix 



Input: a diagram D, m > 2, r > 1, pi, . . . ,Pr £ F^. 
Output: the matrix M associated to C{D; mpi, . . . , mpr). 
Property: computations over ¥p. 

M < — {x"/ : {a, 13) e D}; 

/m < — a one-to-one correspondence from A4 to {1, . . . , if^D}; 
C < — {(fc, dx, dy) : I < k < r,dx + dy < m}; 
fc < — a one-to-one correspondence from C to {1, . . . , }; 
for each {x°'y^, {k,dx,dy)) G A4 x C do 

9 ' — giPk); 

M{fM{x''y^)Jc{k,dx,dy))^g; 
end for each 
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return M; 



Example 16. Let us compute M for D = i)ia£l(3,2, 1), m^'^ = 2^^ pi = (0,0), P2 = (2, 1). The set of 
monomials 

M = {l,x,y,x'^,xy,y'^} 
will be ordered by /m as above, the set of conditions will also be ordered as indicated: 

C = {(1, 0, 0), (1, 1, 0), (1, 0, 1), (2, 0, 0), (2, 1, 0), (2, 0, 1)}. 
Now, for example, take x'^ and (1, 1,0). The polynomial 

d 



9 = 



dx 



(0,0) = 



will be inserted into M in the fj^ix^) = 4th row and the /c(l,l,0) = 2nd column. For the same 
monomial and condition (2, 1,0) we will have 

dx'' 



dx 

inserted into M[4, 5]. The entire matrix is equal to 



(2,1) 



1 








1 











1 





2 


1 











1 


1 





1 











4 


4 














2 


1 


2 











1 





1 



M 



For a matrix M over Fp, let RANk(M) denote the rank of M computed, for example, by using the 
Gauss elimination. 

Algorithm NS 

Input: TO > 2, r > 1, a diagram D, a number of tries t>\. 

Output: NON-SPECIAL or NOT DECIDED, 

NON-SPECIAL implies that L{D;m^^) is non-special. 
Property: computations over Fp. 

repeat t times 

(pi, . . . ,pr) < — randomly chosen points in Fp; 
M < — i-MATRix(£),TO,r,pi, . . . ,Pr); 
k < — RANK(M); 

if fc = min{#Z),r(™^^)} then return NON-SPECIAL; 
end repeat 
return NOT DECIDED; 



Example 17. We will have Ns(2, 2, t)tO0(3, 2, l),t) = not decided, since the system £(0tog(3, 2, 1); 2, 2) = 

£(2; 2, 2) is special. 

Taking D = 0iafl(2, 1), to = 1 and r = 3 we will obtain NS(1, 3, £)ia0(2, l),t) = NON-SPECIAL if and 
only if the algorithm chooses pi,P2,P3 not lying on a line. For a non-special system, the result NON- 
SPECIAL is much more probable if the number t of tries is big. During computations, it appeared that 

"big" in our case means t > 6. 

We will check non-speciality of all systems C{D; m^'^') for D G V, fixed to and all r > 1. 
Algorithm check 

Input: TO > 2, a set of diagrams V, a number of tries t >1. 
Output: a set ^ C f such that for every G G Q, r > 1, 
C{G;m^'^) is non-special. 

for each D gV do 
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if NS(m,r, £),t) = non-SPECIAL then 

if NS(m, r + 1, D,t) — non-SPECIAL then 

g^gu{D}; 

end if 
end if 
end for each 
return Q; 

Observe that the above algorithm is sufficient to check whether all diagrams in V gives non-special 
systems for a fixed multiplicity. However, running it on the set V (of cardinality 17493) from Example [H] 
would consume too much time. Therefore we will reduce all diagrams from D several times (this should 
decrease the number of diagrams) and check whether they are non-special. I yes, we are done due to 
[Dum 09al Theorem 27]. If no, we must deal with diagrams that reduces to special ones. This will be 
explained in more details after presenting the algorithm. 

Algorithm CH 

Input: TO > 2, a set of diagrams T>, 

a number it > of reductions, 

a number w > of reductions performed on reversed diagrams. 

Output: OK or NOT DECIDED, 

OK implies that for every D E T), r > 1, C{D; m^'') is non-special. 

if u > then 

TZ < — red(to, u, V)\ 
TZ i — check(to, TZ, 6); 

V < — REDOUT(m, u, V, TZ); 
end if 

if u > then 

V < — rcv(r'); 

TZ < — red(to,w,2?); 
TZ < — check(to, TZ, 6); 

V < — REDOUT(m, V, V, TZ); 
end if 

TZ < — check(to,2?, 16); 

itTZ = V then return OK else return not decided; 



Example 18. Let us deal with Ch(3,2?, 8,0) for V from Example fTSl 

V = { c)ia0([8] x«), 5ia0([8] "<'), f.ia0([8] ><i°), £)ia0([8] 5, 1), 
Dia0([8] x«, 6, 2), 5ia0([8] 7, 3), Oia0([8] 7, 5)}. 
Each diagram must be 3-reduced 8 times. 



D 


ledf'iD) 




X8-) 


fia0(8,6,2) 


Oia0([8] 


X9-J 


c)ia0(8,8,6,2) 


Oia0([8] 


xlO-j 


ma0(8,8,8,6,2) 


c)ia0([8]x8 


,5,1) 


c)ia0(8,7,5,2) 


c)ia0([8]x8 


,6,2) 


0100(8,8,6,2) 


c)ia0([8]x8 


,7,3) 


c)ia0(8,8,7,3) 


c)ia0([8]x8 


,7,5) 


c)ia0(8,8,7,5) 



We end up with the set TZ containing 6 diagrams (reducing decreased the number of cases). We can 
check that TZ = CHECK(3, TZ, 6), so we are done and the result is OK. An additional advantage lies in the 
size of matrices, since each TO-reduction decreases the number of rows and columns by ('"^ ) . 
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Example 19. The set 2? from Example [5] contains 17493 diagrams. We will run Ch(5, 2?, 3, 0). So we 
must 5-reduce every diagram in V three times, which gives the set TZ of cardinality 6234. All of these 
diagrams appeared to be non-special, so 

7^ = CHECK(5,7^,6) 

and we are done. 

Example 20. We will present the number of diagrams involved in computing Ch(6,X', 14, 13) for 

V = setbignb(6, 51, 8). 

The set T> contains 5472 diagrams. 4617 of them can be 6-reduced 14 times and we obtain the set TZ of 
2991 diagrams. By checking speciality we obtain that 2832 diagrams from TZ are non-special, while the 
rest is probably special. So in V we have 855 not-reducible diagrams together with 250 that reduces to 
special ones. Now we reverse 1105 diagrams, reduce them (all are reducible) 13 times to obtain the set 
with 562 diagrams. Again not all of them are non-special, we are left with 46 diagrams that reduces to 
24 special ones. 

In the next algorithm we deal with systems >C„(a, 5; m^*") for fixed m, n, a and b. Our aim is to 
identify those r, for which the system is special. 

Algorithm finalnba 

Input: TO > 2, 71 > 0, a, 6 > 0. 
Output: the set £ C N such that 

if r ^ £ then £„(a,5;TO^'") is non-special. 

C< — 0; 

if n = then D < — J)ia0([a + 1] ^(^+i)); 

if n > 2 then D < — i)iag([l]^", . . . , [5]^", [6 + l]^('^+i)); 



m. 

repeat 

A < — NS(m, r, £), 16); 
it A~ NOT DECIDED then C < — C U {r}; 
r < — r — 1; 
until A — NON-SPECIAL; 



+ 1; 



repeat 

A < — NS(m, r, D, 16); 
ii A = NOT DECIDED then C < — £ U {r}; 
r < — r + 1] 
until A = NON-SPECIAL; 

return C; 

Example 21. Let us compute FINALNBA(3, 0, 5, 4). We have D = diag{[6]^^), so #D = 30 and, at the 
beginning, r = 5. In the first step we have 

NS(3, 5, D, 16) = NOT DECIDED, 

since in fact £o(5, 4; 3^^) is special. So we take £ = {5} and compute 

NS(3, 6, D, 16) = NON-SPECIAL. 

Since also 

NS(3, 4, D, 16) = NON-SPECIAL, 

we finish with C = {5}. 

The last group of algorithms checks whether a given system is — 1-special, see |Dum 09al Definition 
3]. We begin with auxiliary algorithms take-line and CREMONA (see |Dum 09cl Theorem 3]). We put 

TAKE-LiNE(£((i; mi, ... , rrir)) ~ C{d — 1; mi — 1, TO2 — 1, TO3, . . . , m^), 
CREMONA(£(c?; TOi, . . . , TOr)) — C{d + k; mi + k, TO2 -I- fc, TO3 -|- k, TO4, . . . , to^) 
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for k — d—mi —7712 — 7713. We will also use SOKT{C{d; mi, . . . , rrir)) to sort multiplicities in non-increasing 
order. 

Algorithm spec 

Input: 7M > 2, 71 > 0, a, 6 > 0, r > 1. 
Output: -1-SPECiAL if £„ (a, 6; 777^'') is — 1-special, 

ERROR otherwise. 
Remark: for a system of curves L we define Ld to be the degree, 
Lrrij to be the j-th multiplicity. 

d< — (77 + 1)6 + a; 

TTio < nb + a; 

e < — edim Cn{a, b; m^^); 
for t = 0, . . . ,b do 

L < — £(d - t; 7770, rTi^^*-, (6 -t)^ ("+!)); 

L < — SORt(L); 

repeat 

if Ld - Ljni - L„i2 < then 

L < TAKE-LINE(i); 

else 

if Ld - - i„i2 - L„i3 < then 

L < CREMONA(i); 

end if 
end if 

L < — SORt(L); 

until Ld — L,ni — Lm2 — > 0; 

if edimL > e then return -1-SPECIAL; 
end for 
return error; 

Example 22. Let us compute SPEC(6, 8, 2, 8, 15), so we must consider £8(2, 8; 6^^^). We begin with planar 
system L = £(74; 66, 6^ ^^ 8^) with edimL = -1. For t = take the system £(74; 66, 6^l^ 8^^), sort 
multiplicities to obtain £(74; 66, 8^^, 6^^^). Then use CREMONA four times to obtain £(42; 34, 8, 6^^^). 
Again, use CREMONA to obtain £(36; 28, 6^^''', 2). By the sequence of CREMONA we transform our sys- 
tem into £(8; 2^-"^^). Since edim£(8; 2^^^) — —1, we pass to the case t ^ 1. Now we begin with 
£(73; 66, 6^-^^, 7^^) and, by CREMONA, transform to £(9; 6, 6, 2, 1^^^). Now we use take-line several 
times to produce £(6; 3, 3, 2, 1^^"^). By CREMONA we finish with £(4; 1^^^) of negative expected dimen- 
sion. For t = 2 we begin with £(72; 66, 6^^^) and transform it to £(6; 6, 6). Then, by take-line, we 
obtain £(0; 0) of non-negative expected dimension, so the answer is -1-SPECIAL. 

4. Implementation and results 

All the presented algorithms have been implemented in Free Pascal and can be downloaded from 
[Dum 09b] . They are divided into two kinds, depending on method of working. The first kind simply 
works on given data, the second one prepares batch files with the list of instructions. For example, the 
implementation of red ("red. pas"), of the first kind, performes sequence reductions on given diagrams 
(loaded from the specified file). The algorithm SETPB (of the second kind) prepares the batch file with 
the following instructions (for SETPb(3, 9)) 

tails 3 8,9,10,x,x rt 
basediag 1 1 7 bt 
gluediags inempty bt rt diag 

The above instructions run tails, which produces all admissible £)iag(8, 9, 10)-tails and stores them 
in the file rt; basediag, which prepares Oia0(l, 2, 3, 4, 5, 6, 7) and stores it in bt; gluediags, which glues 
diagrams from inempty (by default, it contains only the empty diagram), bt and rt. 

All algorithms with names beginning with set are of the second kind, together with CH and finalnba. 
The others are of the first kind. 

All algorithms produces log files, where the necessary information is stored. For example, the part of 
log file for multiplicity 2 contains: 
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xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

setbign23 2 2 

result: all systems L_n(a,b) (2"r) are non-special 
for n>=2, a>=0, b>=5, r>=0 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

Itails (all-h-D-admissible tails) 2 3 
tails loaded: 

1 tails loaded, 
tails foimd: 

3 

4 entries used, 2 tails found, 
job finished: 03:11:41:14 

Itails (all-tL-D-admissible tails) 2 4 
tails loaded: 

3 

2 tails loaded, 
tails foimd: 

3 
4 
1 
2 

11 entries used, 5 tails found, 
job finished: 03:11:41:15 

basediag 4 2 
base diagram: 
4,4 

job finished: 03:11:41:15 

tails (admissible tails) 2 

diagram : 

5,5,x 

tails foimd: 

4 

3 

2 

1 

5 

17 entries used, 5 tails found, 
job finished: 03:11:41:15 

gluediags (glue diagrams) 
25 diagrams produced, 
job finished: 03:11:41:15 

check 

multiplicity: 2 

diag(4,4,4) det <> det <> 
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diag 


:4,4,3) 


det <> 





det <> 




diag 


;4,4,2) 


det <> 





det <> 




diag 


;4,4,1) 


det <> 





det <> 




diag 


:4,4,5) 


det <> 





det <> 




diag 


;3,4,4,4) 


det 


<> 





det 


<> 





diag 


;3,4,4,3) 


det 


<> 





det 


<> 





diag 


;3,4,4,2] 


det 


<> 





det 


<> 





diag 


;3,4,4,1) 


det 


<> 





det 


<> 





diag 


;3,4,4,5) 


det 


<> 





det 


<> 





diag 


;4,4,4,4) 


det 


<> 





det 


<> 





diag 


;4,4,4,3) 


det 


<> 





det 


<> 





diag 


;4,4,4,2) 


det 


<> 





det 


<> 





diag 


;4,4,4,1) 


det 


<> 





det 


<> 





diag 


;4,4,4,5) 


det 


<> 





det 


<> 





diag 


;i,4,4,4) 


det 


<> 





det 


<> 





diag 


;i,4,4,3) 


det 


<> 





det 


<> 





diag 


;i,4,4,2) 


det 


<> 





det 


<> 





diag 


;i.4,4,i) 


det 


<> 





det 


<> 





diag 


:i,4,4,5) 


det 


<> 





det 


<> 





diag 


:2,4,4,4) 


det 


<> 





det 


<> 





diag 


:2,4,4,3] 


det 


<> 





det 


<> 





diag 


;2,4,4,2) 


det 


<> 





det 


<> 





diag 


;2,4,4,1) 


det 


<> 





det 


<> 





diag 


:2,4,4,5) 


det 


<> 





det 


<> 






result: positive, 
non-special: 25, special: 
job finished: 03:11:41:16 

The information which set T) is considered, is stored in the preamble. It is then foUowed by names of 
programs together with additional detailed information. The shortlog files contain only preambles and 
names of programs, while inf olog stores only preambles. The f initlog files contain informations on 
running SPEC. Each program informs when it has terminated (day:hour:minute:second). 
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